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Abstract. In this article we initiate the mathematical study of the dynamics 
of a system of nonlinear partial differential equations modelling the motion 
of incompressible, isothermal and conducting modified bipolar fluids in pres- 
ence of magnetic field. We mainly prove the existence of weak solutions to the 
model. We also prove the existence of a trajectory attractor to the translation 
semigroup acting on the trajectories of the set of weak solutions and that of ex- 
ternal forces. Some results concerning the structure of this trajectory attractor 
are also given. The results from this paper may be useful in the investigation 
of some system of PDEs arising from the coupling of incompressible fluids of 
p-structure and the Maxwell equations. 



1. Introduction 
For homogeneous incompressible fluids, the constitutive law satisfies 

T = -tt1 + T(£(u)), 

where T is the Cauchy stress tensor, u is the velocity of the fluid, and ir is the 
undetermined pressure due to the incompressibility condition, 1 is the identity 
tensor. The argument tensor £(u) of the tensor symmetric- valued function T is 
defined through 

£(u) = i(L + L T ), L = Vu, 

where the T superscript denotes the matrix transpose. Magnetohydrodynamics 
(MHD) is a branch of continuum mechanics which studies the motion of conduct- 
ing fluids in the presence of magnetic fields. The system of Partial Differential 
Equations in MHD are basically obtained through the coupling of the dynamical 
equations of the fluids with the Maxwell's equations which is used to take into 
account the effect of the Lorentz force due to the magnetic field (see for example 
[TB]). They play a fundamental role in Astrophysics, Geophysics, Plasma Physics, 
and in many other areas in applied sciences. In many of these, the MHD flow 
exhibits a turbulent behavior which is amongst the very challenging problems in 
nonlinear science. Because of these facts, MHD has been the object of intensive 
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scientific investigation. Due to the folklore fact that the Navier-Stokes is an ac- 
curate model for the motion of incompressible in many practical situation, most 
of scientists have assumed that the fluids are Newtonian whose reduced stress 
tensor T(£(u)) is a linear function of £(u). 

However, there are many conducting materials that cannot be characterized as 
Newtonian fluids. To describe these media one generally has to use (conducting) 
fluid models that allow T to be a nonlinear function of £(u). Fluids in the latter 
class are called Non-Newtonian fluids. We refer for example to the introduction 
of Biskamp's book [12] for some examples of these Non-Newtonian conducting 
fluids. In [33] and [31], Ladyzhenskaya considered a mathematical model for 
nonlinear fluids whose reduced stress tensor T(£(u)) satisfies 

T(5(u)) = 2(e + ^iVur 2 )^, p > 2. (1) 

Since then, this model has been the object of intensive mathematical analysis 
which has generated several important results. We refer to [23], [36] for some 
relevant examples. In [23] the authors emphasized important reasons for consid- 
ering such model. But the fluids or MHD models finding source from ([I]) do not 
really have a meaning in physics as they do not satisfy some principles of contin- 
uum mechanics and thermodynamics. Necas, Novotny and Silhavy [12], Bellout, 
Bloom and Necas [7] have developed the theory of multipolar viscous fluids which 
was based on the work of Necas and Silhavy [UJ . Their theory is compatible with 
the basic principles of thermodynamics such as the Clausius-Duhem inequality 
and the principle of frame indifference, and their results to date indicate that the 
theory of multipolar fluids may lead to a better understanding of hydrodynamic 
turbulence (see for example [ID]). Bipolar fluids whose reduced stress tensor 
T(£(u))is defined by 

T(£(u)) = 2n (e + |£(u)| 2 )^£ (u) - 2 Kl A£(u), (2) 

form a particular class of multipolar fluids. If 1 < p < 2 then the fluids are said 
to be shear thinning, and shear thickening when 2 < p. 

Throughout our analysis, we suppose the existence of a tensor valued function 



snxn 
sym 

that for some p > 1 and for all I, k, i, j = 1, 2, . . . , n, D, E e 



sym ~r ^symi a scalar function E : ffi™^ — > Kq, and constants V\,v 2 such 



iiixn. 
w sym ' 



S (°) =0 > = ^E(O) = 0, T«(D) = ^,E(D), (3) 

dijduVWEijEu >v 1 (l + D)"- 2 |E| 2 , (4) 
^S(D)<z/ 2 (l + |D| 2 r" 2 . (5) 

Here 

Kym = {D£ R nXn : D ij = D jh i, j = 1, 2, . . . , n}. 
In the present work, we will consider a model of MHD arising from the coupling 
of modified bipolar fluids and the Maxwell equations. More precisely, we assume 
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that the reduced stress tensor of the fluids is given by 

T(£(u)) = 2k T(£(u)) - 2«iA£(u). (6) 

Furthermore, let Q be a simply-connected, and bounded domain of lR n (n = 2, 3) 
such that the boundary dQ is of class C°°. This will ensure the existence of a 
normal vector n at each of its point. In this article we are aiming to give some 
mathematical results related to the following non-autonomous partial differential 
equations: 

' fj* - V • T(£(u)) + u • Vu + fiB x curl B + VP = g(x, t) in (0, T] x Q, 

+ ScurlcurlB + /m- VB-/iB- Vu = in (0,T] x Q 
div u = div B = 0, in (0, T] x Q 
u = TijiUjiii = on (0,T) x dQ, 
B n = curl B x n = on (0, T) x dQ, 
u(0) = u , B(0) = B inQ, 

(7) 

where u = (uf, i — 1, . . . , n), B = (Pjj i = 1, . . . , n) and P are unknown functions 
defined on Q x [0,T], representing, respectively, the fluid velocity, the magnetic 
field and the pressure, at each point of Q x [0,T]. The vector n represent the 
normal to dQ. The constants S and \x are positive constants depending on the 
Reynolds numbers of the fluids and the Hartman number. Finally, t^i is defined 
by 

_ dSjjju) 
Tijl ~ d Xl ■ 

The system (j^j) is a modified MHD equation that we obtained from the cou- 
pling of the Maxwell equations and the non-autonomous dynamical equations of 
a modified isothermal incompressible bipolar fluids. The structure of the nonlin- 
earity of problem (JTJ) makes it as interesting as any nonlinear evolution equations 
of mathematical physics such as the basic MHD or the Navier-Stokes equations. 
In addition to well-known tools from Navier-Stokes of MHD, new tools need to 
be elaborated to handle (j7j). 

When B = then (J2J) reduces to the PDEs describing the motion of isothermal 
incompressible nonlinear bipolar fluids which has been thoroughly investigated 
during the last two decades. Several important results were obtained by promi- 
nent mathematical researchers (see, among others, [8], j9], [11], [38], [39J) in this 
direction of research. For p = 2, k% = 0, (JTj) reduces to MHD equations which 
has been the object of intensive mathematical research since the pioneering work 
of Ladyzhenskaya and Solonnikov [35J. We only cite [52], [50], [25], [22] for few 
relevant examples; the reader can consult [27] for a recent and detailed review. 
Assuming that K\ = Samokhin studied the MHD equations arising from the 
coupling of the Ladyzhenskaya model with the Maxwell equations in [UJ, [46] . 
[B], and [15]. In these papers he proved the existence of weak solution of the 
model for p > 1 + Later on Gunzburger and his collaborators generalized 
the settings of Samokhin by taking a fluid with a stress tensor having a more 



TRAJECTORY ATTRACTOR FOR THE MHD OF NON-NEWTONIAN FLUIDS 4 

general p-structure in [28] and [29]. The authors of the later papers analyzed 
the well-posedness and the control of ([7]) still in the case where K\ — and 

In this work we assume that {k\ ^ 0,B ^ 0,p G (1, 1+ ^^]} we are interested 
in the analysis of the long-time dynamics of the weak solutions to ((7|) which is 
very important for the understanding of the global temporal behavior and the 
physical features (such as the turbulence in hydrodynamics) of the model. We 
refer, for instance, to [3J, [33], and [53] for some results in this direction for the case 
of autonomous Navier-Stokes and other autonomous equations of mathematical 
physics. 

In the third section of this paper, we first give an existence result for the weak 
solutions to ([7]) in both two and three dimensional space and p G (1, More 
precisely, by means of a blending of the Galerkin and monotonicity-compactness 
methods we could prove the following 

Theorem 1.1. Let T and E satisfy (u ; B ) £H,jG L 2 (0,T; V*) and 

pe (1, ^gf]- Then there exists at least one pair (u; B) G L°°(0, T; H)nL 2 (0, T; V) 
such that 

and fl7|) holds in the weak sense. 

The notations used in the above statement will be explained clearly in the next 
two sections. 

The long-time behavior of weak solutions to non-autonomous evolutionary 
PDEs ([7]) is analyzed in the last section. The classical and well-known approach 
explained in the monographs [3], [33], and [53] does not apply here since we are 
considering a non- autonomous system. Moreover, when p G (1, ^r§) then we 
do not have a uniqueness result. Therefore, the classical theory of forward at- 
tractor developed in [T5], [3D], [55] and the concept of pullback attractor initially 
elaborated in [201 EH EE] for non-autonomous dissipative dynamical systems does 
not apply here. To overcome the difficulties associated with the possible non- 
uniqueness of solution in the analysis of dynamical system, on can usually use 
three methods. The first one is the method of generalized flow which was initially 
developed in [5]. The second approach is the theory of multivalued semi- flows 
which has already found in the literature in 1948 (see, for instance, [B]). However, 
this method was applied to the analysis of the asymptotic behavior of PDEs for 
the first time in [2] and extended later on in other work ([TJ, [30])- One can also 
cite the paper [13] which extended the notion of pullback attractor to nonau- 
tonomous and stochastic multivalued dynamical systems. The third method is 
the theory of trajectory attractor which was introduced in [TJJ , [37], [32]. Due 
to its simplicity this method has become very popular and used in many works 
but we only mention [17], [18], [19], [31], and [57] for few relevant results. We 
can find a comparison of the two first and the three methods in [13] and [32J, 
respectively. One can also find a recent review on the three methods in [3J. In 
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the present work we use the third approach and we show the existence of the 
trajectory attractor of weak solutions to (J7|) and give some partial results related 
to its structure. More precisely, we obtain the following which will be presented 
in detail in Section 4. 

Theorem 1.2. Let g be translation bounded in Lf oc (M. + ; V*) (that is, it is trans- 
lation compact in Lf oc W (R + ;Y*) ) , then the translation semigroup > 0} act- 
ing on ^s(go) P ossesses a uniform trajectory attractor As(g )- The set -4e( 90 ) ^ s 
bounded in T\ and compact in ® l ° c . It is strictly invariant: 

§tx4.K( S o) = ^s( 90 ), (8) 
for all t > 0. Furthermore, we have 

Ax(g ) = A,(E(go)), ( 9 ) 

where *4^(£(g )) ^ s ^ e uniform (with respect to g G w(£(g )) ) attractor of {/C+, g G 
w(S(^o)}- 

Moreover, for any set B C ^%i go \ which is bounded in J 7 " we aave 

dist c(0jT . D( _ 4 - 5)) (n 0iT § T B, Tl 0iT K z{go) ) -> as t -> oo, (10) 

dist L 2 (0jT . D( _ 4 5 )) (n ,TSTB, n ,TK^( So) ) -> as t -> oo. (11) 

We should mention that even if we drew our inspiration from [T7|, [23], and 
[5TI] the problem we treated here does not fall in the framework of these main 
references. Besides the usual nonlinear terms of the MHD equations it contains 
another nonlinear term of p-structure which exhibits the non-linear relationships 
between the reduced stress and the rate of strain £(u) of the conducting fluids. 
Because of this, the analysis of the behavior of the MHD model (J7|) tends to be 
much more complicated and subtle than that of the Newtonian MHD equations. 
Hence, we have had to invest much effort to prove many important results which 
do not follow from the analysis in the latter works or we could not find in the 
literature. All of our findings rely on these crucial results and they are presented 
in Section 2. 

To the best of our knowledge the present article is the first to deal with (j7|) 
in the setting {k\ ^ 0, B ^ 0,p G (1, 1 + ^^l} - I n ^ ms sense > many topics and 
problems are still open. Some examples of challenges we may address in future 
research are the existence of weak solution for all values of p, the uniqueness 
of such weak solutions. We may also want to study the behavior of the weak 
solution we obtained in this paper as K\ approaches to but fixing p G (1,2); this 
is a very interesting research topic and may lead to an extension of the results 
obtained in jl?], [IB], [II], [IS], [2B] and [29J. These few examples of research 
topics are taken as an analogy of the problems still unsolved in the mathematical 
theory of generalized Non- Newtonian fluids as reported in [9], [25] and [39J. All 
of these questions are very difficult and beyond the scope of this paper, thus we 
will just limit ourselves with giving a suitable mathematical setting for (J7J) and 
partial results related to the dynamics of the weak solutions. However, we hope 
that our work will find its applications elsewhere. 
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The organization of this article is as follows. We introduce the necessary no- 
tations for the mathematical theory of (JTj) in the next section; these notations 
enabled us to rewrite (J7|) into an abstract evolution equations. In the very same 
section we also give very important results that we could not find in the literature. 
Our principal claims rely very much on these new tools. Section 3 is devoted to 
the proof of the existence of weak solutions to ([7]) whose long-time behavior will 
be investigated in the last section of the article. 

Hereafter we should make the convention that the problem ([7]) is characterized by 
its data: 

(1) physical data: Ko, k±, //, Ux, i>2, S,p, n, Q, 

(2) initial conditions uo and Bo, 

(3) and the external force g. 

Unless the contrary is mentioned, any positive constant (C , Ci, Ci, A, etc) depend 
at most on the data of (JTJ). Moreover, they may change form one term to the 
other. 



2. Preliminary: Notations and hypotheses 

We introduce some notations and background following the mathematical the- 
ory of hydrodynamics (see for instance [39]). For any p G [l,oo), L p (<2) and 
W m,p (Q) are the spaces of functions taking values in M n such that each com- 
ponent belongs to the Lebesgue space L P (Q) and the Sobolev spaces W m,p (Q), 
respectively. For p = 2 we use H m (Q) to describe W m,p (Q). The symbols | ■ | and 
(., .) are the L 2 -norm and L 2 -inner product, respectively. The norm of W m ' p (<2) 
is denoted by || • \\ m ,p- As usual, C XD (Q) is the space of infinitely different iable 
functions having compact support contained in Q. The space Wo' m (<2) is the 
closure of C X, (Q) in W p,m (Q). Now we introduce the following spaces 

Vi = {uG Cg°(Q) : divu = 0}, 

Hi = {u e L 2 (Q) : divu = 0, u ■ n = on dQ] , 

Vi = < u G H 2 (Q) : div u = 0, u = = on dQ 

on 



We also set 



V 2 = {B G Cg°(Q) : div B = 0; B • n = on dQ} 



Ho 



the closure of V2 in L 2 (<2), 



Note that 



V 2 = {B G M 1 (Q) : divB = 0; B • n = on dQ] 



Hi = H 2 . 



The spaces Hj, i = 1,2 are equipped with the scalar product and norm induced 
byL 2 (Q). 
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We endow the space Vi with the norm || • ||i generated by the scalar product 

((u - v))l = i g a»» a» t dx > 

where 

£( U ) = ^(L + L T ), L = Vu. 

It is shown in [UJ that this scalar product generates a norm || • ||i which is 
equivalent to the usual H 2 (Q)-norm on Vi. More precisely, there exists two 
positive constants K\ and K 2 depending only on Q such that 

^l||u|| 2 ,2 < ll u l|l < ^2||u|| 2 ,2) (12) 

for any u G Vi. 

On V 2 we define the scalar product 

((u, v)) 2 = (curl u, curl v), 

which coincides with the usual scalar product of H 1 (Q). 
Let 

V = Vi x V 2 , 
H = Hi x H 2 . 

The spaces H and V have the structure of Hilbert spaces when equipped respec- 
tively with the scalar products 

($, vi/) = (u, v) + (B, C), V$ = (u; B), * = (v; C) G H, 

and 

(($, *)) = 2k 1 ((u, v)) a + S{{B, C)) 2 , V$ = (u; B), * = (v; C) G V. 
The norms on H and V are respectively defined by 

|$| = vVl 2 + |B| 2 , V$ = (u; B) G H, 

and 

||$|| = ^KiHuHf + ^||B||1,V$ = (u;B) G V. 
Remark 2.1. Note that the norm ||$|| is equivalent to the norm defined by 

[m] 2 = \Mi + \\B\\l 

We should also mention that a Poincare'-like inequality holds for the space EI and 
V, that is there exists a positive constant A ( dependong only on Q) such that 

A|u| 2 < ||u|| 2 , (13) 

for any u G V. 
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For any Banach space X we denote by X* its dual space and (</>, u) the value 
of G X* on u G X . Identifying HI with its dual, we have the following Gelfand 
chain 

V C H = H* c V*, 

where each space is densely and compactly embedded into the next one. This 
chain of embedding enables us to write 

(u, v) = (u, v), 

for any u G H and v G V. 

Let Ai be a linear operator defined through the relation 

( Am, v) = a (u, v) = 2*1 / d£ f n ) d£ f^ dXi ue^vev, 

Jd oxk 

Note that 

D(Ai) = {u G V : 3f G H C V* for which a(u, v) = (f, v),Vv G V}, 

and A\ = PA 2 , where P is the orthogonal projection defined on L 2 (D) onto 
HI. The operator Ai is self-adjoint and it follows from Rellich's theorem that 
it is compact on Hi. Therefore, there exists a sequence of positive numbers 
{Xi : i = 1,2,3,...} and a family of smooth function {(pi : i = 1,2,3,...} 
satisfying 

Ai<pi = Xi4>i, (14) 

for any i. We can normalize the family {0j : i = 1,2,3,...} so that they will 
form an orthonormal basis of Hi which is orthogonal in Vi- 

We also introduce a linear operator A2 from V2 taking values into (i.e, 
A 2 G C{Y 2 ,Y* 2 ) ) defined by 

(A 2 B,C) = S((B,C)) 2 , 

for any B, C G Y 2 . The operator A 2 is self-adjoint and compact on H 2 , so as 
before we can find a family of increasing positive numbers {//j : % = 1, 2, 3, ... } 
and a family of smooth functions {ipi : i = 1, 2, 3, . . . } such that 

A 2 ipi = mijji, (15) 

holds for any i. It is known from JSU] that the family {ipi : i = 1, 2, 3, . . . } verifies 
the following problem 

curl curl ipi = fiiipi, 
div ipi = 0, 

(V>< • n) aQ = (curl^i • n) dQ = 0, 
for any i. 

Through Ai and A 2 we can define a linear operator from V to V* by setting 

(A$,*) = (AiU,-v) + (A 2 B,C), (16) 
for any $ = (u; B), * = (v; C) G V. 
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To take into account the p-structure of the fluids we introduce a nonlinear 
mapping A p from V into V* by setting 

{Ap®, *) = / T(£(u)) • £(v)dx, 

JQ 

for any \I> = (u; B), \P = (v; C) G V. We state some important properties of A p 
in the following 

Lemma 2.2. Let p G (1, ^gf] and Ze* T and E safe/?/ Tnen, 

(i) £/ie nonlinear operator A p is monotone; that is, 

(Ap^i -Ap$2,®i -$ 2 > > 0, 
/or an?/ $i, $2 £ V. In particular, we have 

(^$,$}>0, (17) 

/or any $ G V. 

(ii) //$ G L 2 (0,T,V) nL°°(0,T,e), then A p <5> G L 2 (0,T,V*). 

To check the results in the above lemma we need to recall the following results 
whose proofs can be found in |25j. 

Lemma 2.3 (Korn's inequalities). Let 1 < p < 00 and let QcM" be of class 
C 1 . Then there exist two positive constants K % = K p (Q),i = 1,2 such that 

Ki\\u\\ hp < \£{u)\*dxy <K*\\u\\ 1>p , (18) 

for any u G Y 1>p . 

Proof of Lemma \2. 2l It is known from (39] that for any p G (1, 00) there exists 
a positive constant u 3 depending only on the physical data such that for all 
D,EeR££: 

(T(D) — T(E)) • (D — E) > 0, (19) 

and 

|T(D)| 2 < 1/3(1 + |D|) 2p " 2 . (20) 
Therefore, we see from ( Fl9]) that for any u, v G Vi 

(A^u - ApV, u - v> = / [T{£{u)) - T(£(v))) ■ [£{u) - £(v)) dx > 0, 

J Q 

which proves the monotonicity of A p . 

By noticing that T(0) = the estimate ffTTj) is a simple consequence of the last 
inequality. 

For any \P = (v; C), $ = (u; B) G V we have 

{Ap*,®) = f T(£(u)) ■ £(v)dx, 
Jq 
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\{M>*)\< 



\T(£(u))\ 2 dx 



\£(w)\ 2 dx 



Q 



This last estimate along with the discrete Holder's inequality and Korn's inequal- 
ity imply that there exists a constant C > such that 



|(^$,*>|<C7 



\T(£(u))\ 2 dx 



0|C||?+|Vv| 2 ]' 



(21) 



As Vi is continuously embedded in Hq(Q), we see from (121]) that 



\{Ap*,9)\<C 



\T(£(u))\ 2 dx 



l^ll. 



Hence 



Let 



<C / \T(£(u))\ 2 dx. 

JQ 



(22) 



X(u) 



\T(£(u))\ 2 dx. 



To deal with x( u ) we w iU distinguish two cases 
Case 1: 1 <p < 2 

We can deduce from (1201) that there exists a positive constant C such that 

|T(D)| 2 < C(l + |D| 2 )(1 + |D| 2 ) P " 2 , (23) 
from which along with p G (1,2] we deduce that 

X(u) <C + C [ \£(u)\ 2 dx. 
Jq 

Owing to Korn's inequality and the continuous embedding of Vi into Hq(Q) we 
infer that 

x(u)<C||u|| 2 . 

Thus, we derive from the last estimate and (122]) the existence of a positive constant 
C such that 

IIA>$llv* <c + c7||$|| 2 , 

which implies that A p § G L 2 (0,T; V*) if $ G L 2 (0,T; V). The proof of item (ii) 
is finished for the case 1 < p < 2. 
Case 1: 2 < p < 2gf 

Throughout this step we set q = 2p — 2. For u G Vi we infer from Gagliardo- 
Nirenberg's inequality that there exists a constant C > such that 



|Vu| 



Q 



vu\ r dx\ <c||u||;|u 



ol,| 1—a 



(24) 
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with 5 < a < 1 and 



r G 

where 



[2, oo) if n = 2, 
[2, s) if n = 3, 



1 _ 1 1 2a 

s n 2 n 

Let q = r G [2, oo) such that ^ = - + \ — ^ with | < a < 1. This is equivalent 
to saying that 

q(2 + n) - 2n 
qa = i . 

The estimate (|24|) implies that 

T rT 

\Vu\\dt<C sup |u(s)|^- a ) 9 / \\u\\l a dt. 

o ' se[o,T] Jo 

Since p G [2, 2 ^rj], then it is not difficult to check that qa < 2 which enables us 
to apply Holder's inequality and infer that 



rp rp OIL 

[ \Vu\ q q dt < C sup |u(s)| (1 - a) " ( / \\u\\ldt) 2 . 
Jo ' se[o,T] \Jo J 

By using Korn's inequality and the embedding of Vi into Hq(Q) we see from 
and (JHD that 



(25) 



X (u(s))ds < C + C \Vu(s)\ q q ds. (26) 
Jo 

Therefore by plugging (T25|) into the last estimate we get that 



/ X(u(s))ds < C + C sup |u(s)| (1 - a)9 ( f \\u\\ldt) 
Jo se[o,T] \Jo J 

Noticing that 

\u(s)\^< [|u(s)| 2 + |B(s)| 2 ]^ : 



(27) 



and 



rp Q a rp ^ a 

J \\u\\ldtY <cH [2^(3)111 + S\\B(s)\\l]dty 



we see that 

rp rp Q2l 

f x(u(s))ds < C + C sup |$(s)| (1 - a) " ( / ||$(s)|| 2 dt) , (28) 
Jo se[o,T] \Jo J 

from which it follows that if $ G L°°(0, T; H) n L 2 (0, T; V) then ^ G L 2 (0, T; V). 
Hence the claim (ii) of the lemma holds true. 

□ 

Before we proceed further, let us state the following 
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Remark 2.4. We could see from the course of the proof of Lemma 12751 that there 
exist positive constants C\ and G% depending only on p, S and Q such that for 
any function $ G L°°(t, t + 1; H) fl L 2 (t, t + 1; V) and t > 

rt+l rt+1 

J ||^||2<Ciy \\Hs)\\ 2 + C 2 , ifpG(l,2], (29) 



and 



<+1 |IA < ^ll2<C'i sup Ms)\^( f T \ms)\\ 2 dt) +C 2 , iipe{2, 



qa 

~* 2n + Q, 



se[o,T] Vo / n + z 

(30) 

where q = 2p — 2 and ga = g( 2+ ^) 2n _ 

As in the case of Navier-Stokes equations we introduce the well-known trilinear 
form b(u, v, w) to deal with the other nonlinear terms of (jZJ). For any u, v, w G 

C °°(Q), we set 



6(u,v,w) = y Ui^-Wjdx 



s 

where summations over repeated indices are enforced. It is also possible to extend 
the definition of the trilinear &(.,.,.) to larger spaces by exploiting the density 
of Cq 3 (Q)( or Vi and V2) in appropriate space. We will do it very often espe- 
cially each time the trilinear form &(., ., .) is continuous. For instance, It is well- 
known (see, among others, [51]) that the trilinear form 6(u, v,w) is continuous 
on H 1 (Q) x H 1 (Q) x M 1 (Q). Moreover, 

6(u,v,v) = 0, (31) 

b(u, v, w) = -6(u, w, v), (32) 

for any u G V 2 , v, w G M l (Q). Since Vi C V 2 , then (J3T]) and (J32]) are also valid 
for any element u in Vi- 
It is also easy to check that 

|6(u,v,w)| <C|u||Vv||w| L o 0(Q) , (33) 

for u G L 2 (Q), v G Hj(Q), w G C °°(Q), and 

|6(u,v,w)| < C|u| l °o(q)|v||Vw|, (34) 

for uGC ^(Q),vGL 2 (Q),wGHi(Q). 
Now we introduce a trilinear form defined on V x V x V by setting 

fio($i, $2, $3) = 6(ui, u 2 , u 3 ) - //6(Bi, B 2 , us) + ^6(ui, B 2 , B 3 ) - ^6(Bi, u 2 , B 3 ), 

(35) 

for any = (u*; B;) G V, i = 1, 2, 3. 

We collect some properties of Bq in the following lemma. 

Lemma 2.5. (i) For any $i,$ 2 G V, there exists a bilinear form £>($i,$ 2 ) 
taking values in V* such that 

Bo($x, $2, $3) = $2), $3), /or $ 3 G V. (36) 



TRAJECTORY ATTRACTOR FOR THE MHD OF NON-NEWTONIAN FLUIDS 13 

(ii) For any $j G V, i — 1, 2, 3 we have 

$ 2j $ 3 ) = 0, (37) 

£ ($i, $ 2 , $3) = $3^2) (38) 

Furthermore, there exists a constant C depending only on the data of (J7J) 
such that 

((£($!, $ 2 ),$ 3 >| <Cmax(l,/x)||0 3 ||[2|$i| ||$ 2 || + ||$i|||$ 2 |], (39) 

for any $j G V, i = 1,2,3. 

(iii) //$ = (u;B) G L 2 (0, T; V) fl L°°(0, T; H) tfien £($,$) G L 2 (0,T;V*). 

Proof. We split the proof into three parts. 
Proof of item (i): 

It is well known that there exists C > such that 

|6(u,v,w)| < C7| |xx| I Iw| [i, 2 | I wj 1 1,2, 

holds for any u, v,w G EI 1 (Q). Hence by using a discrete version of Holder's 
inequality we have 

|£o($i,$2,$ 3 )| <C r max(l,/i)||$ 1 || 1)2 ||$ 2 || li2 ||$ 3 ||i,2, (40) 

for any $j = (u^Bj) G H 1 (Q). We infer from ( l4"0"j) that i3 is continuous on 
V x V x V. Thus, for any given $ 1; $ 2 G V x V there exists an element £>($i, $ 2 ) 
of V* such that Q35) holds for any $ 3 G V*. 
Proof of item (ii): 

We easily derive from (13*11) and (132]) that 

B Q ($i,$ 2 ,$ 3 ) = 0, 

Boi&u $ 2 , $ 3 ) = -Bo($i, $ 3 , $ 2 ), 

for any $j G V, z = 1,2,3. 
Proof of item (iii): 

From dBHU-dBlII and the embedding H 2 (Q) C L°°(Q) we deduce that 

^3 1 1 2,2 (||u 2 ||i|ui| + ||B 2 || 2 |Bi|) 

+ ||B 3 || 2 (||ux|| 2j2 |B 2 | + ||u 2 ||2,2|Bi|) 
By (1121) and discrete Holder's inequality we derive from the last estimate that 

!(£($!, $ 2 ),<I> 3 >| < L7max(l,/i)||0 3 || ||$ 2 || + ||ui|| 2 ,a|B 2 | + | |u 2 | | 2 ,2|Bi| 
Note that 

I | u l| 1 2,2 1 B2 1 + 1 1 "U-2 1 1 2,2 1 Bi I < 1 1 u. 1 1 1 1 "B2 1 + ||Bi|| 2 |u 2 | + ||u||i|Bi| + ||B 2 || 2 |ui| 
which with the discrete Holder inequality yield 

||ui||2,2|B 2 | + ||u 2 ||2,2|Bi| < ||$i|||$ 2 | + H^alH^il- 
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Consequently, 

$ 2 ),$ 3 >| <Cmax(l,/x)||0 3 ||[2|$i| ||$ 2 || + ||$i|||$ 2 |] • (41) 
From ()4ip we obtain 

|B($,$)|v* <3Cmax(l,//)(|$| ||$||), (42) 

which implies that if $ = (u; B) belongs to L 2 (0, T; V)nL°°(0, T; H) then $) 
is an element of L 2 (0, T; V*). □ 

Thanks to all these preliminary consideration, the problem ([7]) can be rewritten 
in the following abstract form 

§ + Ay + A p y + B(y,y)=g, 

y(0) = y„. 1 j 

where y = (u; B) is a solution of and y = (u ; B ). From now on, we will 
work with (|43p . The remaining part of this work is devoted to the analysis of 
(j4"5|) . In the next two sections we will mainly study the existence of its weak 
solutions and their long-time behavior. 



3. Existence of weak solution 

This section is devoted to the investigation of the existence of solutions of (j43p . 
Before we do so, let us define explicitly the concept of solutions that are of interest 
to us. 

Definition 3.1. Let T > and y G H. A weak solution of (g3]) on [0,T] is a 
function y G L°°(0, T; H) n L 2 (0, T; V) such that 

|a 2 (0,T;n 

and 

y(i) - y(r) + ^ ^ly(s)ds + J A p y{s)ds + ^ B(y(s), y{s))ds = J g(s)ds, 

holds as an equality in V* for almost everywhere t, r G [0, T] with t > t. 

If y G L°°(0, T; H) n L 2 (0, T; V) and ^ G L 2 (0, T; V*), then y G (7(0, T; V*) n 
(7(0, T; HLj), where H u denotes the space EI endowed with the weak topology. 
Therefore the initial condition y(0) = yo G H is meaningful. 

The main result of this section is formulated below. 

Theorem 3.2. Let T and £ satisfy ©-(ED, y G H, g G L 2 (0,T; V*) and p G 
(1, ^j]- T/ien t/iere exists at least one solution of (143!) m t/ie sense of Definition 

to" 
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The proof of this statement is based on a blending of Galerkin approximation 
and compactness-monotonicity method. Several steps are needed to achieve the 
target we are aiming for in the current part of the article. 

Step 1: Galerkin Approximation 

For this step we consider the space V m and the operator P m described as 
follows: 

for any positive integer m we set 

V m = Span{* fc = (fa; fa) : k = 1, 2, 3, . . . , m}, (44) 

and P m the orthogonal projection from V*, H onto V m . It is well known (see, for 
instance, [39]) that 

\P m \c(V*;V*) < 1. (45) 
We look for a continuous function y m (t) taking values in V m such that 

^(t) + P m (Ay m (t)+A p y m (t) + l3(y m (t),y m (t)))=P m g, 

y m (0) = y™, 

where y™ is the orthogonal projection of yo with respect to the scalar product of 
H onto V m . The system (j46]) is a system of ordinary differential equations with 
locally continuous coefficients, thus the existence of a continuous function y m (t) 
on a short interval [0, T m ] is ensured by Peano's theorem. The global existence 
will follow from the a priori estimates. As mentioned in the introduction C will 
describe positive constants depending only on the data (not on m) and which 
may change from one term to the next. 

Step 2: Derivation of a priori estimates 

To prove the compactness of our Galerkin solution which will allow us to pass 
to the limit in fT46|) . we need to derive several crucial estimates. 
Multiplying (ggj by y m {t) yields 

'dy™{t) 
dt 

where we have used the fact that 

(P m B(y m (t),y m (t)),y m (t)) = (B(y m (t),y m (t)),y m (t)), 

= ( thanks to (1371).) 

From this we infer that 
l d\y m (t)\ 2 
2 dt 

Owing to ffTTj) we obtain from the last estimate that 



,y m (t) + (P m (Ay m (t) + A p y m (t),y m (t))} = (P m g(t) } y m (t)} } 



\y m (t)\\ 2 + (A p y m (t),y m (t)) = (g(t),y m (t)) . (47) 



2dt 

Using Young's inequality yields 



^|y m Wr + l|y m (t)ll 2 <ll^)Hv*||y m (t)||. 



j t \y m {t)\ 2 + \\y m {t)\\ 2 <\W% 
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Equivalent ly, 

|y m WI 2 + ri|y m (s)ll 2 ^<|yof + f \\g(s)\\Us. 

Jo Jo 

Owing to (145]) we have 



y m (t)\ 2 + / \\y m (s)\\ 2 ds<C(y ,g,T), (48) 



o 



where 



C7(y ,2,T) = |y | 2 + / g(t)dt, 

Jo 



for any m > 0. Since the constant C(yo, g,T) does not depend on m we have 
T = T 

Step 3: Passage to the limit 

Thanks to (I48p and part (ii) of Lemma 12.21 and part (iii) of Lemma 12.51 we can 
derive that belongs to a bounded set of L 2 (0,T; V*). By a diagonal process 
we can find a subsequence of y m which is not relabelled and a function y such 
that 

y m ->■ y weakly-star in L°°(0, T; H), (49) 

weakly in L 2 (0,T; V*), (50) 

y m ->■ y weakly in L 2 (0, T; V). (51) 

Furthermore, by applying the compactness result in [551 Lemma 5.1] we can check 
that 

y m ->■ y strongly in L 2 (0, T; H) . (52) 



Also, owing to (J50|) we see that 

y m (T) -> y(T) weakly in HI. (53) 

Now to complete the proof of the existence of weak solution we have to pass to 
the limit in fl46|) . 

Since A is a continuous linear mapping from L 2 (0, T; V) into L 2 (0, T; V*), then 
we have that 

P m Ay m ->■ ^y weakly in L 2 (0, T; V*). (54) 

Thanks to part (iii) of Lemma 12.51 and ( 148|) P m B(y m , y m ) belongs to a bounded 
set of L 2 (0,T; V*). Taking advantage of (152]) and (15T1). we will show that 

P m B(y m ,y m ) -+ B(y,y) weakly in L 2 (0,T; V*). (55) 

To this end let 

D = {$ = X (t)^ : X(t) e C c °°(0, T) and fc = 1, 2, . . . }, 
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where {^ k ; k = 1, 2, ... } is defined in (j44|) . It is clear that this set is dense in 
L 2 (0,T;V). Owing to [56, Proposition 21.23], the claim ( 15 5|) is achieved if we 
prove that 

' T (B(y m (s),y m (s)) - B(y, y), * k ) X (s)ds -> 0, 



for any $ = %(t) 1 I'fc G D. For this purpose, we rewrite the last identity in the 
following form 



[ (B(y m (s),y m (s)) - B(y, y), 9 h ) X (8)da = h + h- 
Jo 



where 



h = f T (B(y m (s) - y(s), y m (s)), V k )x(s)ds, 
Jo 

I 2 = [ (B(y(s),y(s)-y m (s))^ k } X (s)ds. 
Jo 

For fixed $ and H the mapping T h- >■ J (jS(H, T), ^}x(s)ds is a continuous linear 
functional on L 2 (0,T; V). Hence by invoking (15T1) I2 converges to as m — » 00. 
Next, we easily derive from (133]) that 

< C 



Qf \y m (s)-y(s)\\\y m (s)\\ds 



(B(y m (s)-y(s),y m (s)),* k ) X (s)ds 

\Jo 

x||$ 

which together with Holder's inequality and (j4"81) imply that 



|L°°([0,T]xQ), 



j{B{y m {s) - y(a), y m ( S )), **>x(a)ds < G ( jT |y m (s) - y( S )| 2 ^ 



x ||$| 



L°°([0,T]xQ)- 



Thanks to (1521) the left hand side of this last inequality will converge to as 
m — > 00. Hence we have prove that 1% converges to as m — > 00 which also 
shows that (155]) holds. 

In view of part (ii) of Lemma 12.21 and (148]) . the nonlinear term A p y m is an 
element of a bounded set of L 2 (0, T; V*). Therefore there exists an element O of 
L 2 (0,T; V*) such that 

A p y m -> J) weakly in L 2 (0, T; V*) , (56) 
as m approaches 00. It remains to be shown that 

n = A p y in L 2 (0,T;V*). (57) 
In vertu of ( 15 01) - (1561) . we can pass to the limit in ( 1461) and find that y satisfies 

^ + Ay + n + B(y,y)=g 
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as an equality in L 2 (0, T; V*). Therefore by an integration by parts we see that 

cT r T i 



(n,y(s))ds = J (S( y ( S )), y ( S ))rf S --(|y(T)|^-|y(0)|^), (5? 

where the functional S is defined by 

E($) =g-A$-B(Q,$). 
Since A p is monotone, then for any $ G L 2 (0,T; V) there holds 

<A>y m > y m > > (Ay m , $ ) + (AA y m - 

But from (|46|) we see that 

(Ay m ,y m ) = J o T ^(y m (s)),y m (s)} - ^(|y m (T)| 2 - |y m (0)| 2 ) 

which along with the former estimate and fl58|) yield (after taking the limit as 
m oo) that 

rT 

(Q-Ap<$>(s),y(s)-<f>(s))ds>0. (59) 



At this juncture, we let $ = y - (3$, for any (3 > and * G L 2 (0,T; V). Thus 
after division by (3 and letting — > 0, (loll leads to 

/ (fi-i P y,*>ds>0, V*eL 2 (0,T;V), 
Jo 

from which the sought convergence (|56p - (|57|) follows. Hence thanks to the con- 
vergence (|49|) . (|50|) . ( I5TT) . (|54|) . (155]) and (|57j) we see after passage to the limit 
that there exists at least a weak solution (in sense of Definition 13. ip y to (H3|) . 

4. Trajectory attractor of the solutions of (HH 



In this section we study the behavior of a weak solution y of (1431) for large time. 
We mainly investigate the existence of an attractor a la Chepyzhov and Vishik. 
For this purpose we closely follow the presentation in [17J which is the pioneering 
article dealing with the uniform trajectory attractor of non-autonomous evolu- 
tionary nonlinear partial differential equations. This method does not require 
that a uniqueness of the solution holds. Hereafter, for a Banach space E we set 

rt+i 

L p b (R + ; E) = {u: u(s) G L'l c (R + ; E), sup / u(s)ds < oo}, 

t>o Jt 

and 

rt+l 

INr,?^,.^ = SU P / \Hs)\\ P E ds i 



for any u G L%(M. + ;E), p G [l,oo]. For any r, t > with t > t we define the 
space T^ oc by setting 

Kc = {u : u G Lf oc (R + ; V) n ^(H+j H), ^ G L 2 oc (M + ; V*) } , 
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which is endowed with the metric defined by 



ill 1 i 1 1 2 i i 1 1 2 

l u llj-+ c = II u IIl2 oc(m+ . v) + ||u|| Lj » (R+;H) + 



du I 



dt 



loc' 

Definition 4.1. By Q l ° c we mean the space J-"^ c endowed with the following 



for any u G T x 

Definition 4. 

convergence topology: 

A sequence u m C T^ oc converges to u in 0+ c as m — > oo if 

u m -> u weakly in L 2 (r,t; V), 

u m -> u weakly-star in L°°(r, t; H), 

//1 1 m //l 1 

weak l y in L 2 (r,t;V*), 

for any compact interval [r,t] C M+. 

We define as the space of functions u defined on R + such that 

<9u 



ill i i 1 1 2 i i 1 1 2 

l u ll^ = II u IIl2(m +; v) + ll u llL^(R +; e) + 



dt 



L2(M + ;V*) 



is finite. Note that J 7 " is a Banach space with norm ||u||j™ and any bounded set 
of T\ is precompact in 0^ c . Let g G L 2 oc (IR + ; V*) be an external force which is 
translation bounded, that is 



1 1 1 1 2 ^ 

Definition 4.2. The synbol space £(<?o) is the closure of the set {go{t + ■), t > 0} 
inLg w (R +) V). 

From now on, we denote by Lf ocw (R + ;W*) the space L 2 oc (IR + ;V*) equipped 
with the weak topology. Under some constraints on g the space £(<7o) en j°ys 
some properties which will be important here. We state them below. 

Lemma 4.3. If g G L 2 oc (IR + ; V) is translation bounded then £(<7o) ° compact 
and complete metric space. 

Proof. If go -Lf oc (lR + ; V) is translation bounded the it follows from [T71 Propo- 
sition 6.8] that go is translation compact in Lf ocw (M. + ;Y*), that is the closure of 
{g (t + -),t > 0} in L 2 oc (M+, V*) is compact in L 2 ' OCjW (R + ; V*). Denote this closure 
by %(go)- From [I7J Lemma 6.3], H(go) as a topological subspace of L 2 oc V*) 
is metrizable and the corresponding metric space is complete. Hence 'H(go) can be 
viewed as a compact and complete metric space. But from [T71 Section 8], H(go) 
coincides with £(<7o) up to a homeomorphism. Therefore, £(<7o) as a topological 
subset of L 2 oc (M + ; V*) is a compact and complete metric space. □ 

Remark 4.4. Before we carry on, we should note that the result of the Lemma 
holds true if we consider R in place of R + . 
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Definition 4.5. For a translation bounded function g G L 2 oc (IR + ; V*) we define 
the trajectory space /C+ as the set of weak solutions of fj4"3]) on any compact 
interval / such that the energy inequality 

e at \y(t)\ 2 - e^|y(r)| 2 + [\\y(s)\\ 2 - \\y(s)\ 2 ] < ^( e * - e^), (60) 

holds for almost every r, t G / with t > r + 1 and 

rh 



Q\(g) = sup sup 

Now let 



\fZ\\g(s + t)\\^ds 



he [1,2] *>o \ e Xh 



/Cs( 9 o) = U 9eS ( 3o )/C+ . 
The space /Cs( 90 ) is a topological space equipped with the topology of Q l ° c . We 
recall some crucial properties of /C+ and /Cs( so ). 

Lemma 4.6. (i) For any yo G H and a translation bounded function g G 
L 2 oc (R + ; V*) ; /C+ zs non-empty and the following energy inequality holds 

e A %(t)| 2 - e^|y(r)| 2 + j(V 0|y( S )|| 2 - A|y( S )| 2 ] < 5^(e* - e^), (61) 

/or almost every t, t G M+ mtt t > r + 1 . 
(ii) T/ie se£ zs closed in (B+ c , £(go)) endowed with the usual topology of 
cartesian product space, and /Cs( 90 ) is closed in @^ c . 

Before we proceed to the proof of Lemma 14.61 we recall the following result 
which is taken from [T7J Lemma 9.2]. 

Lemma 4.7. Let f, x G Lj oc (M> + ) and assume that the following inequality holds: 

poo poo poo 

- &(s)f(s)ds + r) f{s)${s)ds< x(s)$(s)ds, 
Jo Jo Jo 

for any $ G C^°(]R_)_; R + ) and for some ?;6R. Then 



f{t)e* - /(r)e^ < J X (s)e" s ds, 
for almost all t, r G M+ t > r. 

Proof of Lemma \4 ■ 6] Let y™ be the projection of yo on the space V m and let 
y m be a sequence of Galerkin solutions of (T4"51) . We have seen from the proof of 
Theorem 13.21 that for any y G H and g G Lf oc {R + ; V*) there exists at least one 
solution to (TJ3]). Now it remains to prove the energy inequality. To do so let us 
consider an element $ G C£°(t, T) with $ > 0. We have 



J\'{t) \\y m {t)\ - \y{t)\\ 2 dt < J\'{t)\y m {t) - y m \ 2 dt, 



(62) 



which implies that v$'|y m | converges to v^lyl strongly in L 2 (r, T). hence, we 
can extract a subsequence from v^ly" 1 ) (not relabelled) such that it converges to 
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\/5P|y| almost everywhere in [r, T\. Since <&'(t) is continuous on a compact set of 
R then it follows from (T4"8"j) that the sequence of numerical functions $'(t)|y m (t)| 2 
is uniformly bounded. Thus we derive from ( |62|) and the Lebesgue Dominated 
Convergence Theorem that 

- j^'( s )\y^( s )\ 2 ds^- j\\s)\y(s)\ 2 ds. (63) 

As a/ $(s) is uniformly bounded on [r, t] , therefore it follows from ( 148]) that 
$(s)y m (s) belongs to a bounded set of L 2 (t, t; V) and 



V / $(i)y m (s) -> v^Iyl^) weakly in L 2 (r,t; V), 
as m — >■ oo. Moreover, 

^(^lly^ll^^hminf^^^lly" 1 ^)!! 2 ^. (64) 

By taking a test function <3> 6 C^°(r,t) with $ > 0, we obtain from fT4T|) that 

~^ /^ /(s)|y? " (s)|2c?S + /^( s )lly T "( s )N 2 ^ + ^^(^)(^y m (s),y m (5)>^ 

<t>(s)(g(s),y m (s))ds. 

Taking into account (TTT|) and using Young's inequality we get from the last esti- 
mate that 

— ^* ^'(^) ly^(^) -H ^* ^(^) I I < ^(^) 1 1 (65) 

Owing to (I6"3"|) and f[64"j) we pass to the limit in (1651) and obtain the following 
energy inequality 



<S>\s)\y(s)\ 2 ds + J <S>(s)\\y(s)\\ 2 ds < J <S>(s)\\g(s)\\ 2 y ,ds, 
which can be rewritten in the following form 

- I\\s)\y m (s)\ 2 ds + Xj^(s)\y m (s)\ 2 ds 

< jf (\\g (s)\\^ - [\\y(s)\\ 2 - X\y(s)\ 2 ]) ds. 
This estimate along with Lemma 14.71 yields 

|y(t)| V - |y(r)|V T + £ (\\y(s)\\ 2 - \y(s)\ 2 ) e Xs ds < J* \\g(s)\\ 2 Y *e Xs ds, (66) 

for almost all r, t G M+ with t > r. To end the proof of item (i) we can use the 
same argument as in [TTJ Corollary 9.4]. 

Item (ii) can be shown by using the same argument as in [TTl Proposition 8.5] 
(see also [TTl Proposition 8.3]). □ 
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We define a semi-group {E> t ,t > 0} acting on £(<?o) and /Cs( 90 ) by 

Stu(-) = u(t + -),t>0. (67) 

Important facts concerning its action on E(go) and /Cs( 90 ) are stated in the fol- 
lowing lemma. 

Lemma 4.8. Assume that g is translation bounded in Z^ oc (IR + ; V*) (i.e., go is 
translation compact in Lf ocw (M. + ; V*) ). Then 

(i) St is continuous on £(<7o) ^ n ^e topology of L^ oc (R + ; V*). 

(ii) We have S f E(p ) C £(<7o) f or an U t > 0. 

(iii) We /iaue 

§4cj;cJ+,V 9 GEy. (68) 

(iv) For any g G £(<7o) we /iave 

^ + C /C+,, (69) 

which implies that 

St/C s(so) C/C s(go) ,t >0. (70) 

Proof. Proofs of (i) and (ii) are straightforward. 

Item (iii) is a direct consequence of fl6lD . Item (iv) is proved as follows: If 
g G £(<7o) and y G then y(i + •),£ > is a solution of fj43l with initial 
condition y(t) and external force g(t + •); that is, Sty G /C§ tff . □ 

To proceed with our investigation we introduce the concept of attractor of our 
interest. The definitions we state below are taken from fTTl. 



Definition 4.9. A set P C Q l ? c is a uniformly (wrt to g G £(g )) attracting set 
for the family {JC g ,g G £(g )} m the topological space 9+ c if for any bounded 
set B of T\ and B C JCwpo) ^ ne se ^ ^* attracts S t B as t — > oo; that is, for any 
neighborhood O(P) in 0+ c there exists t\ > such that StB C O(P) for any 
t > ti. 

Definition 4.10. A set ^4s(g ) is a uniform (wrt to g G £(<7o)) trajectory attractor 
for the family {JC g ,g G S(o )} in the topological space 0+ c if 

(1) ^4s(go) i s compact, 

(2) -4.£(g ) is strictly invariant, that is, S t A^ go ) = A^( go ) for any t G R+, 

(3) and *4e( 90 ) i s a minimal uniformly attracting set for {/C+, g G S(o )}. 

Before we state the main result of this section we still need to introduce the 
concept of w-limit set. 

Definition 4.11. The cu-limit set of £(oo) is defined through 



t>o l h>t Js(9o) 
where [ • ]s( 90 ) designates the closure in E(g ). 
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For any translation bounded function go G L; oc (R + ; V*), Lemma 14.31 implies 
that the metric space £(go) is compact and complete. Hence owing to the item 
(i) of Lemma 14.81 and well-known theorem on continuous semigroup acting on 
compact and complete metric space (see for example [3], [53J ), § t as a continous 
semigroup acting on £(<7 ) nas a global attractor A in £(<?o)- The set A coincides 
with u(H(go)) which implies that 

S t u(E(g ))=u(E(g )), (72) 

for any t > 0. 

From flTT]) we clearly see that a;(£(go)) C T,(g ), so it makes sense to define the 
analog of {JC^,g G H(g )} on u(T>(g ) by taking g in the smaller space u(J2(g )) 
instead of in £(<?o)- Therefore, we can also define the notion of uniform attractor 
for {K+,g G u(E(g ))} with respect to g G u;(£(p ))- 

Now we formulate the main result of this section. 

Theorem 4.12. Let g be translation bounded in Lf oc (R + ; V*) (that is, it is trans- 
lation compact in Lf ocw (R + ;W*)), then the translation semigroup {E>t,t > 0} act- 
ing on A^2( 90 ) possesses a uniform trajectory attractor As( go ). The set -4.£( 90 ) is 
bounded in J 7 " and compact in ® l ° c . It is strictly invariant: 

St«4.s( 90 ) = As(g ), (73) 
for all t > 0. Furthermore, we have 

-4e( 90 ) = A;(E(go)), (74) 

where Aj(s(g )) ^ s the uniform (with respect to g G u(E(go)) ) attractor of {Kg , g G 

As in classical theory of dynamical system of PDEs, to show the existence of 
a trajectory attractor it is crucial to find a uniform attracting set for {§t,i > 0} 
which is bounded in T\ and compact in Q l ° c . Since Q l ° c is a topological space 
endowed with a weak topology defined in Definition 14.11 it is enough to construct 
a bounded set in T\ which will attract S t M for any bounded set B C K^,y For 
the problem (J7|), the following estimate is the main tool for such a construction. 

Lemma 4.13. Assume that go is translation compact in Lf ocw (M> + ; V*) and let 
q = 2p — 2. Then there exist positive constants Ro and Co depending only on the 
data of ([7]) such that 

IIMta < Co\\y\\l°o (OA . m) e- M + Rv,Vt > 1, ifp G (1,2], (75) 



and 



if P e(2, 



|§*y||^ < C \\y\\ q Loo{om e'^ + i? ,W > 1, (76) 



2ra+6 l 
n+2 \ 



Proof of Lemma We distinguish the case p G (1, 2] and p G (2, ^±51 

Case I: p g (1,2] 



> n+2 
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The proof will be split in several steps. Throughout, y will denote an element of 

£(£(«))) and 9 e ^(9o) ■ 

Estimate for §<y in Lg°(R + ;E[) 

Owing to the Poincare'-like inequality (TT3|) we derive from f l60|) that 



e M \y(t)\ 2 -e A -|y(r)| 2 < 



A 



ie xt -e 



(77) 



holds for almost every r,t 6 I with t > r + 1. Recall that for a fucntion go which 
is translation bounded in Lf oc (R + ;Y*), the following holds (see [T7J Section 8]) 



\9\\h( 



< 



I 1 12 

Ifl'ollx? 



which combined with (I77|) yields 



2 -At 



|y(t + s)| 2 < |y(s-t)| 2 e 
for t > 1 and s > 2. Therefore, we obtain that 



5^(l- e -*), 
A 



|S*y| 



< 



11 1 2 -At , ^A(ffo) 

|j| lL°°(0,l;H) e "T 



A 



(7f 



for any t > 1. 

Estimate for S t y in L^(R+;V) 



Hereafter let us set R\ 



. By integrating ( |6T|) between t and t + 1 and 



invoking (IT8"j) . we obtain that 
/•t+l 

A y |y(s)| 2 rf s < X\\y\\l^ m + Rie xt (e x - 1), 
which along with (I6T|) and (J78|) implies that 
/' + (lly(s)|| 2 - A|y( S )| 2 ) e Xs ds < XR l e x \e x - 1) + i? ie At + ||y 

By vertu of fT79|) and flHOj) we have 
/■t+i 

y ||y(s)||V s ds < 2R 1 (e x - l)e At + R x e xt + (1 + A)||y 



|2 

lz,°°(o,i;: 



(79) 



10) 



L°°(0,1;H)- 



For s G [t,t + 1], we can easily check that e Ai ||y(s)|| 2 < e As | |y(s) 1 1 2 . Hence, we 
infer from flHTl) the existence of positive constants R2 and C2 such that that 



|S*y| 



< R 2 + C 2 e A *IIy||ioo( ,i; 



(82) 



for any t > 1. 

Estimate for S t (gr) in Lg(R + ;V*) 
From (143 1), it is clear that 



<9y 

at 



< C7 (H^lyll 2 , + ||^l p y|| 2 , + ||i3(y,y)|| 2 , + ||g || 2 ,) 
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With the definition of A (see (USD), the estimate (J2SD of Remark EH and (@2J) we 
derive from fl83l) that 



rt+l 

Jt 



dy(s) 



dt 



rt+l rt+l 

ds<C sup |y(s)| 2 / \\y{s)\\ 2 ds + C \\g (s)\\ 2 ds 

s£[t,t+i\ Jt Jt 



+C 



t+i 



|y( S )|| 2 + C7 snp \y(s)\ 2 + C 
se[t,t+i] 



(84) 

where the constants involved in this estimate depend only on the data of (JZJ). We 
infer from this last estimate, (ITS"]) and fl82j) that there exists positive constants C 3 
and i?3 such that 



■t+i 



<9y(s) 



9* 



< ^3||y||x,«j(o,l;H) + -^3 



55) 



for any t > 1. 

Now, the proof of the lemma follows from (jTSjl . (1821 and fl85l) if p G (1,2]. 
Case II: p g (2, fg] 



We should first not that the estimates (j78|) and (jSZj) hold true for the case p < 2 
and p > 2, the main difference lies in th estimating of §j (%)■ Throughout this 
step I = p — 1 which will belong to (1, 2) provided that p G (2, 2 ^M^]. Because of 
this fact, we can derive from Holder's inequality, ( 17S|) and (I8"2"j) that 



qAt 



4( K +> 



< i? 4 + C 4 e"~||y 



qAt 



< flg + Cse™ y 



lL°°(0,l;H)' 



lL°°(0,l;H)' 



for any t > 1. Owing to these last inequalities we infer from f l5U|) . (I7S|) and 
that 



l^y(s)||v^s < ^6 + C 6 ||y| 



2 

Lg°(0,T;H)l 



9 19 
I 2 



qAt 



ILP°(0,T;] 



which easily implies that 

r-t+1 



HAy(s)llv*^<i?7 + C 5 ||y| 



With these estimates we can derive from ( EH ) that 

2 



L?°(0,T;H) e 2 ' 



t+1 



<9y(s) 



at 



ds 



q\t 

< C 8 e -V 



|yllL°°(Q,l;H) + "^8 



Now the proof of (jTbj) follows from all of these estimates. 



□ 



Proof of Theorem \4-l%\ First of all, Lemma fl~3l implies that the symbol space 
£(<7o) is a compact and complete metric space. Secondly, we see from Lemma 
14.61 that the family {1C+,g G £(g )} is (0+ c , £(g ))-closed. Moreover, it follows 
from Lemma [4.81 that the translation semigroup § t is continuous on T,(g ), and 
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Sf£(<?o) C S(po) f° r au ^ > 0, and the family {JC+,g G £(#o)} is translation 
coordinated (that is, it verifies (169]) ). Thanks to Lemma [4.131 the ball Br 

B Ro = {ueJ*:\\u\\ n <2R }, 

is a uniformly attracting set for {JC+,g G £(#0)} in @+ c - The set I?# is compact 
in 0+ c and it is obviously bounded in T\. Thus the conditions of [UJ Theorems 
3.1] can be applied to our case, and this proves our theorem. □ 

In the next assertion we will study the structure of the uniform trajectory 
attractor of Theorem 14.121 To this end we should extend all the concepts we 
have introduced throughout this section to the case where the external force go 
has an extension to R. This is possible thanks to the invariance property (I72]l 
and arguing exactly as in [T7J Section 4]. Then, we denote by J-/ oc and J 7 " the 
analog of J-^ c and J 7 " by taking R in place of R + . In a similar way we also define 
Qioc p Qr a f unc ti on u defined on R we designate by Il + u its restriction to R + . 

Definition 4.14. The complete symbol space Z(g ) is the set of all functions 
C G Lf oc (R; V*) such that U+C(t + ■) the restriction of ((t + ■) to R+ belongs to 

w(£(<7o)) f° r an Y t G M. 

For any £ G Z(go), the kernel is the set of all weak solutions y to fj4*3]) with 
external force £ such that y G Jioc and y satisfies the energy inequality fl60|) for 
almost everywhere t, r G R with t > r + 1. For a translation bounded function 
2oGLL(R;V*) we set 

K z(so) = |J Kf. 
Cez(go) 

As in [T7] we also have the following result. 



Theorem 4.15. Under the assumptions of Theorem 4-12 , 

Ax{ go ) = A>(E( So ))> = n + U K ^ = U + K z(go)- ( 86 ) 

<ez(g ) 

The kernel is not empty for any ( G Z(g ), and the set K Z ( go )) is bounded 
in J 70- and compact in ® l ° c . 

Proof. We omit the proof of this result as it is very similar to the proof of [17] 
Theorem 4.1]. □ 

From the construction we made in Section [2] that A is self-adjoint and thanks 
to Rellich's theorem it is compact on HI. Therefore we can define the fractional 
power operators A 5 with domain D(A S ) for 5 G R. For any reflexive Banach 
space X and Y let 

u(s), s G [0, T] : u G L p ° (0, T;X),^e L Pl (0, T; Y) 

which is a Banach space when endowed with the graph norm. Taking 5 G [0, 1/2) 
implies that the embedding EI C D(A~ S ) is compact and D(A~ S ) C V* is continu- 
ous. Therefore, owing to [5U Corollary 4] the space W^, 2(0, T, H, V*) is compactly 
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embedded in C(0, T; D(A s )). Since V C D(A 5 ) is compact then by invoking 
[361 Lemma 5.1] we see that W 2 , 2 (0, T, V, V*) is compact in L 2 (0,T; D(A 5 )). The 
consequence of these facts and Theorem 14.151 are given below. 

Corollary 4.16. For any set B C ^s(go) w ^ c ^ ^ s bounded in T\ we have 

dist c . (0iT;D( ^-« )) (rio,TSTB, Yl 0tT K z{go) ) as t ->■ oo, (87) 
dist L 2 (0jT;D( _ 45)) (n 0iT § T B, n , T IK z(9o) ) as t ->■ oo. (88) 

In the corollary ITo/ru stands for the restriction to [0, T] of a function u defined 
on R (or R + ) and dist(X, F) denotes the distance between the set X and Y. 
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